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ABSTRACT 


High  latitude  propagation  phenomenon  such  as  scatter  from 
the  aurora  and  sporadic-E  propagation  can  alter  the  communication 
characteristics  of  the  high  latitude  meteor  communication  chan¬ 
nel.  This  report  describes  the  results  of  several  experiments 
conducted  for  the  purpose  of  characterizing  the  fading  and 
multipath  profile  of  the  high  latitude  meteor  communication 
channel.  The  effect  of  high  latitude  propagation  phenomena  on 
the  capacity  and  reliability  of  meteor  burst  communications  is 


discussed . 


EXECUTIVE  SUMMARY 


Meteor  burst  communications  has  several  potential  advan¬ 
tages  over  conventional  HF  communications  at  high  latitudes. 
Because  it  does  not  depend  upon  ionospheric  reflection,  it  is  not 
subject  to  the  diurnal,  seasonal  and  solar  induced  changes  in  the 
maximum  usable  frequency  (MUF)  as  is  HF.  Network  control  is 
therefore  much  simpler  since  single  frequency  operation  can  be 
used  and  no  adaptive  frequency  management  is  required.  Meteor 
burst,  due  to  its  small  reflecting  area,  has  a  certain  degree  of 
inherent  AJ  and  LPI  to  ground  based  interceptors  and  jammers. 
Meteor  burst  operates  at  higher  frequencies  than  HF  (40  to 
100  MHz  for  meteor  burst  versus  2  to  30  for  HF)  making  it  much 
less  vulnerable  to  ionospheric  absorption  events.  For  these 
reasons ,.  meteor  burst  is  being  considered  or  is  being  applied  to 
several  high  latitude  communication  missions. 

There  are  several  high  latitude  propagation  effects  which 
can  alter  the  communication  characteristics  of  the  channel;  some 
of  which  can  be  considered  beneficial,  while  some  could  be  detri¬ 
mental.  The  objective  of  the  program  described  in  this  report 
was  to  design  a  relatively  inexpensive  experiment  to  characterize 
the  various  propagation  mechanisms  in  addition  to  meteors  which 
are  observed  at  high  latitudes.  The  primary  objective  was  to 
determine  the  Doppler  and  multipath  spreads  of  the  various 
propagation  modes,  and  to  determine  their  potential  effect  on 
conventional  meteor  communication.  In  particular,  three  propa¬ 
gation  mechanisms  which  result  in  multipath  and  fading  at  VHF  at 
high  latitudes  have  been  characterized:  one  is  multipath  due  to 
multiple  and/or  warped  meteor  trails,  another  is  sporadic-E 
scatter,  and  the  third  is  scatter  from  the  aurora. 


i)  The  simultaneous  occurrence  of  more  than  one  meteor  on  a 
single  link  can,  for  a  short  period  of  time,  produce  multipath 
spreads  as  high  as  1000  us;  however,  for  a  typical  link  configur¬ 
ation,  only  on  1  to  2  percent  of  all  meteor  trails  will  there  be 
significant  multipath  interference  due  to  multiple  trails.  On 
occasion,  individual  trails  can  warp  so  that  multiple  reflection 
points  are  produced  creating  potential  multipath  up  to  a  maximum 
of  about  5  us.  This  occurrence  is  fairly  rare  and  even  then  the 
coherence  bandwidth  limit  due  to  multipath  is  several  hundred 
kb/s.  Warping  of  individual  trails  due  to  high  altitude  wind 
shears  with  multipath  spreads  less  than  1  us  is  a  primary  source 
of  fading  observed  on  many  meteor  trails.  This  and  several  other 
studies  have  shown  that  the  meteor  channel  is  primarily  signal 
strength  limited  as  opposed  to  multipath  limited  channels,  such 
as  the  troposcatter  channel .  Advanced  adaptive  techniques  which 
adapt  throughput  to  channel  signal  strength,  such  as  the  variable 
data  rate  modem,  can  operate  on  the  meteor  channel  and  can 
greatly  increase  throughput. 

ii)  Sporadic-E  layer  scatter  at  low  VHF  frequencies  (30  to 
60  MHz)  occurs  frequently  at  high  latitudes  resulting  in  a 
continuous  strong  signal  which  exhibits  little  multipath 
dispersion.  Due  to  the  unreliability  and  low  predictability  of 
the  sporadic-E  channel,  as  opposed  to  the  intermittent  yet  more 
reliable  meteor  channel,  this  channel  is  best  considered  as  a 
"bonus"  channel  rather  than  as  the  primary  survivable,  reliable 
backup  channel.  There  are  several  ramifications  of  the  frequent 
occurrence  of  sporadic-E  at  high  latitudes.  Protocols  for  meteor 
communication  systems  must  be  designed  to  operate  in  a  high 
contention  environment  when  a  continuous  channel  to  many  links 
may  exist  simultaneously  for  long  periods  of  time.  Simple  meteor 
burst  networks  operate  on  the  principle  that  the  duty  cycle  on 
any  link  is  so  low  that  collisions  between  transmissions  are 


-xi- 


rare.  The  sporadic-E  channel  due  to  its  large  reflection  area, 
has  much  lower  jam  and  intercept  resistance  than  does  the  meteor 
channel,  subjecting  the  meteor  network  to  potential  interference 
from  unintentional  or  intentional  signal  sources  operating  on 
adjacent  frequencies.  The  primary  multipath  source  on  the 
sporadic-E  channel  is  ironically  the  occurrence  of  meteors  off  to 
the  side  of  the  great  circle  path  producing  intermittent  multi- 
path  spreads  up  to  1000  us  on  longer  links.  The  Doppler  spread 
of  sporadic-E  channels  is  on  the  order  of  10  to  30  Hz  or  less. 
Figures  ES-1  and  ES-2  show  the  multipath  and  Doppler  profile  of  a 
typical  sporadic-E  channel.  We  observe  the  occurrence  of  meteors 
in  Figure  ES-1  in  frames  21  ,  33,  3  8,  44  and  49. 

iii)  Scatter  from  the  aurora  is  a  localized  phenomenon  occurring 
when  the  aurora  and  the  meteor  link  are  oriented  such  that  condi¬ 
tions  for  specular  scatter  from  elongated  irregularities  in  the 
ionosphere  are  satisfied.  When  present,  auroral  scatter  can 
represent  an  additional  background  interference  source,  or  on 
occasion,  depending  on  the  link  power,  antenna  patterns  and 
density  of  the  auroral  ionization,  it  can  create  a  high  multipath 
(multipath  spreads  up  to  1000  us  or  more),  fast  fading  (200  Hz 
Doppler  spread),  high  signal  strength,  continuous  channel. 
Auroral  scatter  is  much  more  localized  than  is  sporadic-E  since 
the  requirements  for  field  aligned  auroral  scatter,  outlined  in 
the  companion  report  by  Malaga  [1986],  are  much  more  restrictive 
than  for  simple  layered  sporadic-E. 


Figure  ES-3  shows  a  series  of  multipath  snapshots  of  the 
Anchorage-Bethel,  Alaska  meteor  channel  acquired  every  six 
seconds  on  September  1  9,  1985.  In  Figure  ES-3  the  large  peaks  in 
frames  47,  53,  and  54  represent  large  meteor  echoes  and  the 
remainder  is  due  to  auroral  scatter.  The  five-minute  average  rms 
multipath  spread  in  Figure  ES-3  is  801  us.  Figure  ES-4  shows  the 
Doppler  spread.  In  Figure  ES-4,  the  component  centered  at  zero 


RELATIVE 

AMPLITUDE 


Figure  ES-1  Series  of  multipath  snapshots  showing  Sporadic-E  event  on  1  Jul^ 

1985  at  0300  UT.  Snapshots  were  acquired  every  6  seconds.  Multi- 
path  due  to  meteors  is  clearly  visible  in  snapshots  21,  33,  38,  44 
and  49. 
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Fiqure  ES-4  Composite  Doppler  profile  for  period  shown  in  Figure  ES-3. 

Component  centered  on  zero  shift  corresponds  to  meteor  propagation 
and  the  shifted  component  is  due  to  the  auroral  scatter. 


frequency  is  due  to  scatter  for  meteors  and  the  component  cen¬ 
tered  at  about  80  Hz  is  due  to  scatter  from  the  aurora.  These 
figures  show  how  auroral  scatter  can  transform  the  benign  meteor 
scatter  channel  into  a  fast  fading  high  multipath  channel. 

The  degree,  if  any,  to  which  auroral  scatter  will  disrupt  a 
given  meteor  communication  link  depends  upon  the  antenna  pat¬ 
terns,  transmitter  power,  and  the  location,  orientation  and  ioni¬ 
zation  density  of  the  aurora  relative  to  the  link.  In  the  exper¬ 
iment,  we  have  demonstrated  that  the  potential  exists  for  disrup¬ 
tion.  Communication  using  coherent  modulation  at  conventional 
meteor  burst  data  rates  (4  to  8  kbps)  in  the  fast  fading,  high 
multipath  environment  shown  in  Figures  ES-3  and  ES-4  would  be 
virtually  impossible. 

At  a  minimum,  protocols  for  a  network  operating  in  auroral 
regions  must  be  designed  to  allow  operation  of  a  network  even  if 
there  are  continuous,  yet  potentially  unusable  channels  contend¬ 
ing  for  the  network.  To  actually  communicate  over  the  disturbed 
auroral  channel  two  techniques  can  be  used.  Advanced  equaliza¬ 
tion  techniques  using  implicit  diversity  (such  as  those  applied 
to  the  mul tipath- limited  troposcatter  channel)  can  be  used  to 
exploit  the  characteristics  of  the  multipath  spread.  A  simpler 
technique  is  to  reduce  the  data  rate  to  less  than  the  multipath 
spread  and  use  coding  with  non-coherent  modulation  to  maintain 
comrnun  ication . 

Taking  into  account  the  undesirable  occasional  disruption 
due  to  auroral  scatter,  meteor  burst  still  has  many  advantages 
over  HF  which  can  be  even  more  effected  by  these  same  phenomena. 
Auroral  scatter  degradation  at  meteor  burst  frequencies,  may  be 
absorption  blackout  conditions  at  HF . 

The  work  described  in  this  and  the  companion  theoretical 
development  by  Malaga  represent  an  important  first  step  in  scien¬ 
tifically  evaluating  the  reliability  of  meteor  burst  communica- 
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tion  at  high  latitudes.  The  research  results  must  now  be  applied 
to  practical  communication  systems.  A  new  meteor  burst  comunica- 
tion  network  will  become  operational  in  1986  to  provide  austere 
backup  communication  between  remote  radar  sites  in  Alaska  and  the 
regional  control  center  in  Anchorage.  A  natural  and  important 
follow-on  to  the  experiment  described  in  this  report  should  be  to 
collect  data  from  the  Alaskan  Air  Command  on  the  performance  of 
the  various  links  in  the  network  to  assess  which  links  are  domi¬ 
nated  by  sporadic-E  channels,  and  which  links,  if  any,  are 
limited  by  auroral  scatter  effects.  Data  should  be  compared  to 
the  predictions  of  Malaga  to  calibrate  and  validate  the  model. 

Links  showing  auroral  scatter  performance  degradation 
should  be  instrumented  to  determine  the  multipath  and  Doppler 
spread.  Using  the  validated  model  and  data  from  the  experiment, 
techniques  to  mitigate  the  degradation  could  be  developed 
including  adaptive  antenna  steering,  reduced  data  rate  operation, 
increased  operating  frequency,  and  adaptive  message  routing. 
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SECTION  1 
INTRODUCTION 


Meteor  burst  commun  ication  uses  the  ionized  trails  of  mete¬ 
ors  for  long  distance  digital  communication.  At  high  latitudes 
meteor  burst  has  several  potential  advantages  over  HF  for  long 
haul  digital  communication  appl  ica  t  ions .  Because  meteor  burst 
operates  at  higher  frequencies,  40  to  100  MHz  versus  2  to  30  MHz 
for  HF,  it  is  much  less  sensitive  to  polar  cap  absorption  events 
which  can  create  HF  blackout  conditions.  Meteor  burst  can  oper¬ 
ate  at  a  constant  frequency  unlike  HF  which  requires  frequency 
management  due  to  the  diurnal  and  seasonal  changes  in  the  maximum 
usable  frequency  (MUF).  Meteor  burst  is,  however,  not  totally 
immune  to  high  latitude  effects  and  in  order  to  accurately  assess 
the  reliability  of  meteor  burst  as  an  alternative  to  HF  at  high 
latitudes,  a  thorough  understanding  of  the  characteristics  of  the 
high  latitude  meteor  channel  is  required. 


Previous  high  latitude  meteor  burst  efforts  [Maynard,  1972; 
Ostergaard,  1985;  Cannon,  1985]  have  focused  on  determining  the 
statistics  of  meteor  arrivals,  durations,  waiting  time  distri¬ 
butions  and  the  effect  of  polar  cap  absorption  events  on  these 
statistics.  The  fading  and  multipath  profile  of  the  high  lati¬ 
tude  channel  has  been  assumed  to  be  similar  to  the  mid-latitude 
channel.  Recent  experiments  indicate  that  this  assumption  is  not 
totally  valid  because  it  neglects  the  effects  of  frequently 
occurring  high  latitude  propagation  phenomena  such  as  scatter 
from  the  aurora  and  sporadic-E  layer  propagation.  This  report 
discusses  the  results  of  several  experiments  in  which  the  fading 
and  multipath  profiles  of  the  high  latitude  meteor  channel  were 
characterized.  The  results  are  combined  with  the  knowledge  of 
the  mid-latitude  channel  to  develop  a  latitude  dependent  multi- 
path  and  Doppler  profile  for  the  meteor  channel. 
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Development  of  a  multipath  model  for  the  high  latitude 
meteor  channel  begins  in  Section  2  with  a  review  of  the  charac¬ 
terization  of  the  middle  latitude  meteor  channel.  Multipath  on 
the  middle  latitude  meteor  channel  occurs  primarily  due  to  the 
occasional  occurrence  of  more  than  one  meteor  or  due  to  multiple 
reflections  from  a  single  trail  which  has  warped  due  to  high 
altitude  wind  shears.  Results  from  several  middle  latitude 
multipath  characterization  experiments  are  presented. 

Section  3  of  this  report  discusses  several  propagation 
mechanisms  occurring  primarily  at  high  latitudes  which  can  de¬ 
grade  or  enhance  meteor  communication.  These  mechanisms  incude 
scatter  from  the  aurora,  and  sporadic-E  related  phenomena.  The 
theory  of  these  propagation  mechanisms  is  presented  in  Malaga 
[1986].  Two  experiments  to  characterize  multipath  phenomena  of 
VHF  channels  at  high  latitudes  are  described. 

The  first  experiment,  described  in  Section  4,  used  the 
large  bandwidth  (330  KHz)  allocation  on  the  1260  km  USAF  high 
latitude  meteor  burst  test  bed  from  Sondrestram  AB,  Greenland  to 
Thule  AB ,  Greenland.  High  resolution  multipath  measurements  were 
performed  on  sporadic-E  and  meteor  channels  in  the  45  MHz  band. 

The  second  experiment  described  in  Section  5,  used  a  dedi¬ 
cated  link  (640  km)  from  Anchorage,  Alaska,  to  Bethel,  Alaska, 
operating  at  42.4  MHz  to  perform  lower  resolution  multipath 

measurements  on  sporadic-E,  auroral,  and  meteor  channels.  The 

✓ 

Alaska  link  represented  the  primary  research  tool  in  the  multi- 
path  measurement  program. 

Finally,  Section  6  discusses  the  ramifications  of  the 
findings  on  the  reliability  of  meteor  burst  communication  at  high 
latitudes . 


SECTION  2 

THE  MULTIPATH  STRUCTURE  OF  THE  METEOR  CHANNEL 


2.1  DIFFUSE  MULTIPATH  FROM  A  SINGLE  METEOR  TRAIL 

Meteor  channels  form  when  the  ionized  trails  of  meteors 
entering  the  atmosphere  satisfy  the  following  geometric  require¬ 
ments  shown  in  Figure  2-1. 


1.  The  point  of  tangency  of  a  line  defining  the  center  of 

a  cylindrical  meteor  trail  to  an  ellipsoid  of  revolu¬ 
tion  with  foci  at  the  receiver  and  transmitter  lies 
within  the  volume  common  to  both  the  receive  and 

transmit  antenna  patterns. 

2.  The  region  of  the  cylinder  defining  the  meteor  trail 

which  contains  ionization  fills  at  least  one-half  of 
the  Principal  Fresnel  Zone.  The  Principal  Fresnel 

Zone  is  the  region  of  the  trail  such  that  the  distance 
from  transmitter  to  receiver  through  a  point  on  the 
trail  differs  by  less  than  one-half  wavelength  from 
the  distance  to  the  point  of  tangency.  All  scatterers 
within  this  region  contribute  constructively  to  the 

received  signal  while  scatterers  within  higher  order 

Fresnel  zones,  which  are  much  smaller,  contribute 

alternately  constructively  and  destructively  to  the 

received  signal.  The  length  of  the  Fresnel  Zone  is 

dependent  on  link  distance  and  frequency.  Typical 
meteor  trails  are  15  to  50  km  in  length  and  the 

Principal  Fresnel  Zone  is  about  1  to  4  km  at  50  MHz. 


At  the  time  of  formation,  meteor  trails  can  be  modeled  as 
lines  or  narrow  cylinders  of  electrons.  Immediately  after  forma¬ 
tion,  the  trails  begin  to  expand  due  to  diffusion.  In  most 
trails  the  electron  density  is  low  enough  so  that  each  electron 
is  assumed  to  scatter  independently  of  all  other  electrons  and 
the  received  power  is  computed  by  integrating  the  contribution 
from  all  scatterers  over  the  entire  trail.  This  type  of  trail  is 
referred  to  as  an  underdense  trail.  Closed  form  solutions  for 
the  power  scattered  by  underdense  trails  and  its  decay  as  the 
trail  expands  have  been  derived  by  Eshlemman  [1955]. 


Trail  Showing  Primary  and  Secondary  Fresnel 


Underdense  trails  form  rapidly  and  then  decay  exponentially 
with  time  constants  on  the  order  of  200  to  400  ms  at  50  MHz 
[Sugar,  1964].  From  a  multipath  viewpoint,  the  underdense  meteor 
channel  is  very  benign,  limited  only  by  the  several  Fresnel  zone 
diameter  of  the  trail  as  it  expands.  Each  Fresnel  zone  into 
which  the  trail  radially  expands  adds  multipath  corresponding  to 
one-half  wavelength  path  delay  or  about  10  ns  at  50  Mhz.  It  is 
the  integration  of  the  Gaussian  radial  electron  density  of  the 
trail  as  it  diffuses  radially  into  higher  order  Fresnel  zones 
that  causes  the  exponential  decay  of  the  total  scattered  power. 

The  coherence  bandwidth,  defined  as  the  inverse  of  the 
delay  spread,  is  a  measure  of  the  maximum  data  rate  that  can  be 
transmitted  without  intersymbol  interference.  If  symbols  are 
transmitted  at  a  rate  greater  than  the  coherence  bandwidth,  total 
overlap  of  adjacent  symbols  will  occur.  For  practical  systems 
which  do  not  employ  equalization,  the  maximum  symbol  rate  to 
avoid  intersymbol  interference  should  be  about  one-tenth  the 
coherence  bandwidth.  The  coherence  bandwidth  of  underdense 
meteor  channels  is  on  the  order  of  several  megahertz  allowing 
symbol  rates  without  intersymbol  interference  as  high  as 
1  million  symbols/sec.  Unfortunately,  the  several  megahertz 
coherence  bandwidth  of  the  underdense  trails  is  difficult  to 
exploit  since  in  most  meteor  communication  applications,  the 
channel  is  signal  power  limited  (due  to  the  150  to  160  dB  path 
loss)  rather  than  multipath  limited.  Many  trails  can  support 
data  rates  in  the  hundreds  of  kilobits  per  second  for  a  short 
time  if  fully  adaptive  communication  systems  are  used  [Weitzen, 
1983,  Weitzen,  et  al . ,  1984]. 

As  the  trail  density  and  received  signal  power  increases, 
the  model  for  the  trails  begins  to  change  as  the  assumption  that 
all  electrons  scatter  independently  begins  to  break  down.  The 
model  for  the  received  power  becomes  complicated  as  the  electrons 
begin  to  interact  and  only  rough  closed  form  approximations  to 


the  received  power  exist  [Hines  and  Forsythe,  1957],  These 
trails  are  referred  to  as  overdense.  As  the  trails  become 
denser,  the  lifetime  above  a  given  signal  threshold  tends  to 
increase  so  that  high  altitude  wind  shears  have  time  to  effect 
the  trail  geometry.  These  winds  can  act  on  an  initially  linear 
meteor  trail  shown  in  Figure  2-1  until  it  becomes  warped  and 
distorted  as  shown  in  Figure  2-2.  Fresnel  zones  are  superimposed 
on  both  Figures  2-1  and  2-2  for  reference.  Note  in  Figure  2-2 
that  portions  of  the  trail  which  were  initially  in  the  Principal 
Fresnel  zone  contributing  constructively  to  the  received  signal 
are  now  in  higher  order  Fresnel  zones  contributing  destructively. 
The  result  of  the  time  changing  trail  distortion  is  fading,  al¬ 
though  the  delay  spread  of  the  trail  is  still  very  small  because 
the  majority  of  the  trail  lies  within  several  Fresnel  zones. 

2.2  DISCRETE  MULTIPATH  FROM  A  DISTORTED  METEOR  TRAIL 

On  occasion,  a  trail  can  distort  so  that  several  regions  in 
the  trail  meet  the  requirements  for  specular  reflections  outlined 
earlier.  This  phenomenon,  which  occurs  much  more  frequently  for 
backscatter  geometries  than  forward  scatter,  was  first  described 
by  Manning  [1959].  Sugar,  et  al.  [1960]  extended  the  backscatter 
work  of  Manning  to  forward  scatter  and  estimated  some  worst  case 
multipath  spreads.  As  an  example,  consider  a  1200  km  link  with  a 
20  km  meteor  trail  oriented  normal  to  the  great  circle  path  with 
one  end  of  the  trail  at  the  mid-point  of  the  path.  Assume  that 
a  second  reflection  point  occurs  at  the  far  end  of  the  trail. 
Using  simple  geometry,  the  worst  case  delay  between  the  two  com¬ 
ponents  would  be  3.7  us.  Most  spreads  are  far  less  than  the 
worst  case.  These  spreads  coincide  with  the  experimental  meas¬ 
urements  of  Weitzen  [1983]  and  Carpenter  and  Ochs  [1962]. 
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2.3  DISCRETE  MULTIPATH  FROM  MULTIPLE  METEOR  TRAILS 


An  event  that  could  cause  more  severe  multipath  inter¬ 
ference  is  the  simultaneous  occurrence  of  two  or  more  meteors 
within  the  common  volume  of  the  receiver- transmitter  antenna 
pattern.  Relative  delays  beween  components  of  1  ms  are  possible. 

We  now  attempt  to  predict  how  often  this  event  might  occur. 
Let  Tl  be  the  lifetime  of  a  trail  defined  as  the  interval  during 
which  the  power  scattered  by  the  trail  exceeds  an  arbitrary 
signal  threshold.  A  multiple  trail  event  occurs  when,  during  the 
lifetime  of  a  given  trail,  Tl,  one  or  more  additional  trails 
occur  which  exceed  the  received  signal  threshold. 

Meteor  trails  arrive  randomly  so  that  the  probability  of 
exactly  k  trails  arriving  within  a  time  T,  Pk(t),  can  be  de¬ 
scribed  by  a  Poisson  process  as 

-LCT 

Pk(T)  =  ^-yr~  acT)k  (2.1) 

where  Lc  is  the  average  arrival  rate  ( trails/sec ) of  meteor  trails 
whose  scattered  signal  strength  exceeds  the  minimum  signal 
strength  required  for  reliable  communications.  To  aid  the  analy¬ 
sis,  we  will  assume  that  all  meteor  trails  decay  according  to  the 
underdense  model.  This  assumption  is  often  used  because  it 
allows  the  memoryless  property  of  the  exponential  decay  to  be 
used  to  develop  closed  form  solutions.  Using  this  assumption, 
the  distribution  of  trail  lifetimes,  Tl»  was  theorized  by  Sugar 
[1964]  to  be  exponentially  distributed  as 


P(Tl>T)  =  e  T  (2.2) 
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D  is  the  average  diffusion  constant  =  8  nr/sec 


The  probability  that  during  the  lifetime  TL  of  a  trail,  one 
or  more  trails  occur  which  exceed  the  communication  threshold  is 
given  by 


Pfone  or  more  trails  trail  with  lifetime  )  = 


-  (LcTJ  ~LcTL 

L  ITT  e  “  1  ~ 


-LcTL 


(2.3) 
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It  can  be  seen  from  Equation  (2.3)  that  for  trails  with  short 
duration,  the  probability  of  multiple  trail  occurrence  is  very 


small.  The  probability  of  multipath  interference,  P 


m 
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determined  by  integrating  the  density  conditioned  on  lifetime 


Equation  (2.3)  over 
Equation  (2.2)  as 


the  exponential  density  of  lifetimes 
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As  an  example,  assume  that  the  arrival  rate  Lc  of  meteor  trails 
which  exceed  the  communication  threshold  is  2.5  trails  per  minute 
(0.0416  trails/sec).  Let  the  average  lifetime  t  of  the  trails  be 
400  ms.  From  Equation  (2.4),  there  is  a  1.6  percent  probability 
that  there  will  be  some  degradation  in  the  error  probability  due 
to  multipath  interference. 

2.4  REVIEW  OF  PREVIOUS  MULTIPATH  EXPERIMENTS 

Several  experimental  efforts  have  been  undertaken  to  char¬ 
acterize  the  multipath  profile  of  the  meteor  channel  at  mid¬ 
latitudes.  Carpenter  and  Ochs  [1962],  conducted  the  most  exten¬ 
sive  experimental  program  using  an  800  kW  (standard  meteor 
communication  systems  generally  have  transmitters  with  less  than 
1  kW  output)  peak  pulsed  power  forward  scatter  radar  with  a  3  us 
pulse  duration  operating  at  a  frequency  of  41  MHz.  With  the  high 
pulsed  power  signal,  a  weak  component  due  to  ionospheric  scatter 
was  observed  in  addition  to  meteor  propagation.  The  primary 
multipath  mechanism  observed  by  Carpenter  and  Ochs  was  the 
simultaneous  occurrence  of  more  than  one  meteor  within  the  common 
volume.  With  the  800  kW  transmitter,  the  average  time  between 
the  arrival  of  meteors  was  on  the  order  of  several  hundred 
milliseconds.  They  observed  that  on  about  6  percent  of  the 
meteor  channels,  a  second  discrete  component  less  than  6  dB  down 
from  the  main  component  was  present.  Multipath  spreads  varied 
from  the  resolution  of  the  equipment  up  to  the  maximum  1  ms 
resolution. 

In  order  to  predict  the  effect  of  two-trail  multipath  on 
actual  communication  systems,  they  interpolated  the  results  of 
the  800  kW  experiment  to  a  1  kW  communication  link  and  compared 
the  predictions  to  the  performance  of  an  actual  system  [Sugar, 
et  al.,  I960].  Depending  on  the  value  of  several  constants,  they 
predicted  that  multipath  would  occur  from  1  percent  to  3  percent 
of  the  time.  The  predictions  compare  well  to  their  experimental 
measurements  and  to  Equation  2.4. 
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On  many  occasions,  they  observed  enhanced  propagation  due 
to  sporadic-E  or  enhanced  E-layer  propagation.  The  majority  of 
these  events  were  characterized  by  one  discrete  component;  how¬ 
ever,  during  one  of  these  events,  two  discrete  layers  with  20  us 
relative  delay  were  observed.  In  addition  to  multipath  due  to 
discrete  layers,  the  combination  of  sporadic-E  located  near  the 
mid-point  of  the  great  circle  path  and  large  meteors  several 
hundred  kilometers  to  the  side  of  the  mid-point  path  caused 
multipath  with  spreads  up  to  the  900  us  resolution  of  the  experi¬ 
ment  . 

On  one  occasion  during  a  solar  disturbance,  scattering  from 
what  is  believed  to  be  the  aurora  was  observed.  Continuous 
spreads  of  at  least  500  us  were  observed  for  two  hours. 

Weitzen,  et  al.  [1984]  performed  an  experiment  using  tele¬ 
vision  signals  at  50  MHz  (Channel  2)  as  a  wideband  probing 
waveform,  (6  MHz  bandwidth)  extending  the  work  of  Grossi  and 
Javed  [1971].  One  microsecond  resolution  was  obtained  and  little 
multipath  was  observed,  confirming  the  findings  of  Carpenter  and 
Ochs . 


SECTION  3 

EXPERIMENTS  TO  CHARACTERIZE  THE  HIGH-LATITUDE  VHF  CHANNEL 
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The  characteristics  of  the  VHF  meteor  channel  at  high  lati¬ 
tudes  appear,  from  preliminary  results  obtained  using  the  USAF 
High  Latitude  Meteor  Test  Bed  [Ostergaard,  et  al.,  1985],  to  be 
similar  to  the  charac teri sties  at  middle  latitudes.  However 
other  propagation  mechanisms,  such  as  scatter  from  the  aurora  and 
sporadic-E  propagation,  occur  much  more  frequently  at  high  lati¬ 
tudes  than  at  middle  latitudes,  creating  the  potential  for  ex¬ 
tended  periods  when  meteors  are  not  the  dominant  propagation  mode 
on  a  VHF  1  ink  . 

We  have  taken  a  two-step  approach  to  characterizing  the 
effects  of  these  mechanisms  on  meteor  communication  links.  Prior 
to  an  experimental  program  Malaga  [1986  ]  developed  the  theory  of 
VHF  scatter  from  the  aurora  so  we  would  know  what  to  expect  from 
the  experiment.  Next,  a  two  phase  experimental  program  was 
developed  which  is  described  in  the  next  sections. 


3.1  EXPERIMENTAL  GOALS  AND  REQUIREMENTS 

The  objective  of  the  experimental  program  was  to  develop 
and  implement  a  relatively  inexpensive  technique  to  characterize 
multipath  and  Doppler  spreads  on  a  number  of  different  types  of 
channels.  To  characterize  the  channel  completely,  multipath 
spread,  Doppler  spread,  and  signal  strength  statistics  are  re¬ 
quired.  In  the  experiment  we  decided  signal  strength  statistics, 
though  desirable,  would  be  secondary  to  the  primary  multipath  and 
Doppler  spread  statistics.  This  decision  was  based  on  the  added 
costs  required  to  accurately  calibrate  and  maintain  calibration 
on  the  receiver,  the  potential  inability  to  maintain  constant 
transmitter  power  on  the  Alaska  multipath  link,  and  the  fact  that 
detailed  signal  strength  statistics  for  meteor  scatter  propaga- 
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tion  and  sporadic-E  propagation  are  currently  being  collected  on 
the  high  latitude  testbed.  A  further  concern  during  the  experi¬ 
ment  design  was  the  requirement  that  the  equipment  be  totally 
automated  to  a  point  where  the  only  service  required  would  be,  at 
most,  bi-weekly  maintenance  to  change  data  tapes. 

A  real-time  RAKE  receiver  or  a  very  high  power  wide-band 
pulsed  radar  are  two  techniques  which  are  best  suited  to  multi- 
path  probing.  Both  techniques  are,  however,  expensive  and  diffi¬ 
cult  to  implement:  the  real-time  RAKE,  due  to  the  complicated 
hardware  required  to  perform  correlations  in  real-time,  and  the 
pulsed  radar  due  to  the  very  high  power  (Carpenter  and  Ochs  used 
800  kW)  transmitter  required.  Acquiring  data  and  processing  it 
at  a  later  time  is  a  third  technique  which  has  the  advantage  of 
being  much  less  expensive  than  the  real-time  RAKE  receiver.  The 
disadvantages  of  this  technique  are  large  data  storage  require¬ 
ments,  potentially  high  data  analysis  costs,  and  inability  to 
provide  on-site  monitoring  of  multipath  and  Doppler  spreads. 

The' compromise  solution  selected  for  the  experiment  was  to 
use  the  signal  processing  capability  of  the  Data  Precision  D-6000 
signal  analyzer  to  preprocess  the  data  on-site.  Data  would  be 
collected  during  a  portion  of  each  hour  and  processed  during  the 
remaining  time  providing  near  real-time  monitoring  of  channel 
conditions.  Only  preprocessed  data  would  be  retained  for  analy¬ 
sis  at  a  later  time.  The  D-6000  analyzers  were  available  due  to 
their  extensive  use  in  the  high  latitude  program  as  were  the 
cassette  tape  drives  and  the  HP-85  controller.  While  other 
equipment  might  have  been  more  optimal  for  the  application,  the 
cost  savings  incurred  using  the  existing  equipment  and  the  exper¬ 
ience  of  RADC/EEPS  personnel  in  the  design  of  experiments  based 
on  this  equipment  weighed  heavily  in  the  decision. 

A  preliminary  experiment  was  designed  to  test  the  correla¬ 
tion,  data  processing,  and  control  algorithms  on  the  existing 
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USAF  High  Latitude  Meteor  Test  Bed  in  Greenland.  A  secondary 
objective  was  to  use  the  330  KHz  bandwidth  available  on  the  link 
to  investigate  the  fine-structure  of  sporadic-E  multipath  and/or 
multiple  meteor  trail  occurrences.  The  results  of  this  experi¬ 
ment  are  described  in  Section  4.  Figure  3-1  shows  the  1260 
Greenland  link  from  Sondrestrom  to  Thule. 

With  the  lessons  of  the  preliminary  experiment,  a  large 
scale  experiment  was  undertaken  to  characterize  the  multipath  and 
Doppler  spreads  of  auroral  scatter,  sporadic-E  and  meteor  chan¬ 
nels.  Cost  constraints  required  scaling  down  the  experiment  to 
use  leased  transmission  time  from  a  meteor  burst  common  carrier 
company  in  Anchorage  instead  of  a  dedicated  link  for  transmission 
services.  The  results,  nevertheless,  were  successful  and  are 
described  in  Section  5.  Figure  3-2  shows  the  Primary  Link  from 
Anchorage  to  Bethel,  Alaska. 


Location  of  Auroral  Scatter  Multipath  Measurement  Link  in  Alaska 


SECTION  4 

PRELIMINARY  MULTIPATH  EXPERIMENT  IN  GREENLAND 


A  preliminary  experiment  to  develop  and  test  channel 
probing  algorithms  and  waveforms  and  to  characterize  the 
sporadic-E  channel  with  a  330  KHz  waveform  was  performed  in 
Greenland  using  the  U.S.  Air  Force  High  Latitude  Meteor  Burst 
test  bed  with  transmitter  at  Sondrestrom  Air  Base  and  receiver  at 
Thule  Air  Base.  This  link  is  one  of  the  few  locations  where  wide 
transmission  bandwidths  are  still  available  in  the  low  VHF  band 
(crowding  in  the  low  VHF  band  would  prevent  an  experiment  such  as 
the  Carpenter  and  Ochs  experiment  from  being  repeated  in  the  con¬ 
tinental  United  States  in  the  foreseeable  future).  The  link  is 
well  above  the  auroral  oval  and  auroral  scatter  is  not  generally 
observed  on  the  link;  however,  high  latitude  sporadic-E  propaga¬ 
tion  is  frequently  observed. 

4.1  EQUIPMENT  FOR  GREENLAND  EXPERIMENT 

Transmit  equipment  located  at  Sondrestrom  Air  Base, 
Greenland,  shown  in  Figure  4-la,  consisted  of  an  800  Watt  synthe¬ 
sized  transmitter  with  330  kHz  bandwidth  allocations  at  45  MHz, 
65  MHz  and  104  MHz.  The  frequency  allocation  at  45  MHz  was 
selected  to  maximize  the  observation  of  sporadic-E  and  meteors. 
Five  element  Yagi  antennas,  1.5  wavelengths  above  the  ground  were 
used  at  both  the  transmit  and  receive  locations.  The  peak  gain 
of  the  antennas  was  about  10.5  dBi . 

A  maximal  length  cyclic  pseudorandom  (PR)  sequence  was 
transmitted  as  the  probing  waveform  (instead  of  a  pulsed  waveform 
as  in  the  Carpenter  and  Ochs  experiment),  to  provide  processing 
gain  to  offset  the  lower  transmitter  power.  The  tradeoff  in 
using  the  PR  sequence  is  that  more  complicated  receiver  process¬ 
ing  is  required.  The  length  of  the  PR  sequence  was  chosen  so  as 
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Figure  4-la  Configuration  of  Transmit  Figure  4-lb  Receive  Equipment  Configu- 

Equipment  for  Preliminary  ration  for  Greenland  Multi- 

Multipath  Experiment  in  path  Experiment. 

Green!  arid  . 


to  optimize  the  processing  gain  and  the  maximum  multipath  spread 
that  can  be  resolved  while  minimizing  the  degradation  in 
processing  gain  due  to  frequency  offsets  between  transmitter  and 
receiver  and  Doppler  shifts  when  noncoherent  detection  is  used. 
Equation  (4.1)  relates  the  processing  gain  to  the  net  frequency 
offset,  sequence  length  and  chip  rate. 


sin  it  vLT  2 

Processing  Gain  =  L( — — — ) 


where 

L  =  Sequence  length 

v  =  Frequency  offset  between  transmitter  and  receiver 
including  Doppler  shifts  and  offsets 

Tc  =  Chip  period 

Ec  =  Chip  energy 

Ng  =  Noise  density 

Table  4-1  lists  the  characteristics  of  the  probing  waveform 
selected  for  the  experiment.  Assuming  calibrated  frequency  syn¬ 
thesizers  at  receiver  and  transmitter  so  that  the  worst  case 
transmitter  -  receiver  frequency  offset  was  50  Hz,  and  assuming 
at  most  100  Hz  Doppler  shift  on  meteor  trails  and  sporadic-E,  the 
worst  case  processing  gain  degradation,  assuming  high  signal  to 
noise  ratio,  is  calculated  from  Equation  (4.1)  as 

Gain  degradation  =  -0.318  dB 


2Nn  itvLT  ; 

1  +  _ 2.( _ 5 _ ) 

LE  vs in  (  ttvLT  )  j 
c  c 


(4.1) 
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Table  4-1  Waveform 
Greenl and 

Characteristics  for 
Exper imen  t 

Sequence  Length  (L) 

127  Bits 

Maximum  Frequency  Offset 
Between  Rx  and  Tx 

±50  Hz 

Maximum  Doppler 

Shift  Expected 

±100  Hz 

Chip  Rate 

128  kb/s 

Chip  Period  (Tc)  and 

Multipath  Resolution 

7.8125  ps 

Sequence  Period  or  992.29  ps 

Maximum  Multipath  Resolution 

Processing  Gain 
(Ideal ) 

21  dB 

Worst  Case  Processing  Gain 
Degradation  (high  SNR) 

-0.318  dB 

Number  of  Samples  per 
Sequence  Chip 

8 

Sampling  Rate 

1.024  x  106  samples/s 

Bi  t  s/Sampl  e 

9 
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The  receive  equipment  located  at  Thule  Air  Base,  Greenland 
(see  Figure  4-lb)  consisted  of  custom  designed  low  noise  receiver 
with  10.7  MHz  IF,  quadrature  IF  to  baseband  downconver ter ,  data 
acquisition  system,  data  cassette  mass  storage  device  and  HP-85 
controller.  The  receiver  had  two  outputs,  an  IF  output  which  was 
non-coherently  quadrature  downconverted  to  baseband,  and  a  log- 
amplitude  detector  calibrated  to  measure  received  signal  power  in 
the  IF  bandwidth  of  330  kHz.  A  strip  chart  recorder  was  connec¬ 
ted  to  the  log-ampl  itude  detector  to  record  long  term  events  and 
to  aid  in  the  determination  of  the  multipath  mechanism.  RF  and 
digital  portions  of  the  receive  equipment  were  isolated  in  sepa¬ 
rate  rooms  to  minimize  interference.  The  data  acquisition  system 
consisted  of  a  Data  Precision  D-6000  digital  oscilloscope/signal 
analyzer  with  an  18  MHz  dual  channel  A/D  converter. 

Data  was  acquired  whenever  the  output  of  the  log-amplitude 
detector,  used  as  a  trigger,  exceeded  a  given  threshold.  Quadra¬ 
ture  data  was  9  bit  quantized  at  a  rate  of  8  samples  per  PR 
sequence  chip  (1.024  Msamples/sec)  and  stored  in  D-6000  memory. 
In  order  to  observe  possible  changes  in  the  multipath  profile 
during  the  lifetime  of  a  meteor  trail,  each  time  the  analyzer  was 
triggered,  a  second  trigger  was  forced  about  850  ms  later.  After 
acquisition  of  two  quadrature  data  records,  post-processing  was 
performed.  Processing  involved  quadrature  correlation  of  the 
data  samples  with  a  stored  version  of  the  transmitted  PR  sequence 
followed  by  computation  of  the  square  of  the  I  and  Q  correlation 
components.  The  two  processed  waveforms  and  the  quadrature 
correlations  were  then  stored  on  tape  for  future  analysis. 
During  post-processing  which  required  about  four  minutes  for  the 
two  complex  waveforms,  data  acquisition  was  disabled.  Real-time 
matched  filter  correlation  would  have  been  a  more  efficient 
technique  allowing  for  triggering  on  smaller  meteors  (the  20  dB 
processing  gain  is  not  achieved  until  the  post  trigger  proces¬ 
sing)  and  reducing  the  time  required  for  processing;  however,  the 
cost  of  the  equipment  would  have  been  prohibitive. 


4.2  RESULTS  FROM  THE  GREENLAND  MULTIPATH  EXPERIMENT 


Data  was  collected  for  approximately  10  days  during  March 
1985,  after  which  interference  with  sensitive  air  traffic  control 
equipment  at  Sondrestrom  AB  by  harmonics  of  the  wideband  trans¬ 
mitted  signal  forced  early  termination  of  the  experiment.  Data 
from  the  experiment  was  preprocessed  according  to  the  procedure 
described  earlier  and  stored  on  tape  cassette.  The  final  step  in 
the  data  analysis  which  occurred  after  the  experiment  was  con¬ 
cluded  involved  measuring  and  determining  statistics  of  the 
multipath.  See  auroral  scatter  experiment  for  details  of  the 
process.  Figure  4-2  shows  an  example  of  the  correlations  in 

which  one  discrete  meteoric  scatterer  was  observed.  Figure  4-3 
shows  an  example  in  which  two  discrete  correlation  peaks  were 
observed.  Using  the  strip  chart  recording  as  a  guide,  it  was 
determined  that  the  mechanism  causing  the  multipath  is  most 
likely  two  meteor  trails.  The  hundred  microsecond  delay  between 
components  eliminated  trail  warping  as  a  possible  mechanism  and 
the  strip  chart  showed  no  sign  of  background  ionospheric  or 
sporadic-E  scatter.  Figure  4-4  shows  two  very  closely  spaced 

components  (less  than  8  us  apart)  which  are  most  likely  two 
closely  spaced  meteors  but  could  possibly  be  due  to  warping  of  a 
single  trail.  Figure  4-5  shows  the  single  discrete  echo  observed 
during  a  sporadic-E  event.  Figure  4-6a  shows  the  probable  occur¬ 
rence  of  a  meteor  trail  and  sporadic-E.  In  Figure  4-6b,  taken 
850  ms  later,  the  trail  has  decayed  away. 

Due  to  the  long  time  required  to  process  the  two  1016  (127 
chip  sequence  x  8  samples  per  chip)  point  complex  correlation 
waveforms,  sequences  of  snapshots  showing  the  time  history  of 

sporadic-E  events  were  not  possible.  This  defficiency  was 

corrected  in  the  Alaska  experiment.  Due  to  early  termination  of 
the  experiment,  inadequate  data  were  obtained  for  detailed  sta¬ 
tistical  analysis;  however,  several  qualitative  conclusions  from 
the  limited  data  could  be  drawn. 


Figure  4-4  Two  Closely  Spaced  Components  Most  Likely  Caused  By  Two  Meteors  But 
Possibly  Due  to  Warping  of  a  Single  Trail.  Spacing  of  the  peaks  is 
1  ess  than  8  gsec  . 
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About  1500  pairs  of  waveforms  (trigger  followed  by  second 
trigger  850  ms  later)  were  obtained  during  the  10  day  period. 

.  About  550  of  the  pairs  were  discarded  for  various  reasons,  in¬ 

cluding  false  triggers  and  times  during  which  the  transmitter  was 
not  operational . 

Of  the  remaining  926  records,  885  were  due  to  meteors  and 
^  the  remaining  41  triggers  occurred  during  sporadic-E  events.  Of 

the  885  meteor  records,  approximately  2  percent  showed  discrete 
multipath  due  to  multiple  meteors  or  trail  warping  within  the  8 
to  992  us  resolution  of  the  experiment.  The  average  spacing 
between  discrete  components  was  226  ys.  Sporadic-E  was  observed 
on  several  occasions  and  meteors  combined  with  sporadic-E  were 
observed.  Table  4-2  summarizes  the  Greenland  experiments. 

While  the  data  from  this  experiment  was  somewhat  inconclu¬ 
sive  due  to  the  early  termination  of  the  experiment,  the  many 
lessons  learned  during  the  experiment  were  applied  to  the  Alaska 
experiment  and  contributed  to  its  success.  These  lessons  were: 

1.  Collect  data  during  events,  process  later  rather 
than  on  the  fly.  This  will  capture  data  during 
infrequent  events. 

2.  Add  a  CW  probe  period  to  perform  FFT  measurements 
to  observe  fading  statistics. 

3.  Use  less  samples  per  chip  to  reduce  processing 
,  t ime . 

4.  Collect  more  data  records  during  a  shorter  period 
of  time  to  provide  time  history  of  events.  Use 
snapshots  every  few  seconds  rather  than  trigger¬ 
ing  on  signals. 

5.  Compute  histograms  of  received  signal  level. 

Eliminate  the  second  trigger  850  ms  after 
first.  It  was  of  very  little  use  and  would  have 
been  of  use  only  if  we  had  greater  resolution 
(less  than  1  us )  to  observe  trail  warping. 
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Table  4-2  Summary  Greenland  Multipath  Experiment 


Total  Pairs  Acquired 

1482 

Pairs  Discarded 

556 

Valid  Data  Pairs 

926 

Number  of  Meteor  Records 

885 

Number  of  Meteor  Records  Showing 

Multipath  (meteors) 

17 

(1.92  %) 

Average  Spacing  between 

Discrepancy  Components 

226  us 

Number  of  Sporadic-E  Records 

41 

Number  of  Sporadic-E  Records 
Showing  Multipath 

4 

(9.7%  ) 
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SECTION  5 

ALASKA  MULTIPATH  EXPERIMENT 


The  primary  experiment  was  conducted  in  Alaska  near  the 
auroral  region  with  the  objective  to  observe  multipath  effects 
due  to  the  combination  of  meteors,  sporadic-E  and  scatter  from 
the  aurora.  The  maximum  frequency  allocation  that  could  be 
obtained  in  Alaska  was  15  kHz,  using  an  existing  meteor  burst 
common  carrier  company.  A  high  power  transmitter  (10  kW)  with 
16  dBi  transmitter  antenna  was  used  to  increase  the  probability 
of  observing  faint  background  ionospheric  and  auroral  scatter.  A 
ten  minute  period  during  each  hour  was  leased  from  the  common 
carrier  service  during  which  a  probing  waveform  was  transmitted. 
The  probing  sequence  consisted  of  four  minutes  of  CW  for  Doppler 
measurements  and  RSL  histograms,  five  minutes  during  which  a 
cyclic  redundant  31  bit  PR  sequence  was  transmitted  at  8  kbps  for 
multipath  measurements,  and  a  one  minute  silent  period  for  noise 
measurements.  The  transmitter  located  near  Anchorage  operated  on 
an  assigned  carrier  frequency  of  42.4  MHz  in  a  15  kHz  allocation. 

5.1  EQUIPMENT  FOR  ALASKA  EXPERIMENT 

A  receive  site  was  required  in  which  both  auroral  scatter 
and  sporadic-E,  in  addition  to  meteors,  were  observable  using  the 
existing  common  carrier  meteor  burst  transmission  facility  in 
Anchorage.  Based  on  several  prediction  models  [Malaga,  1986]  in 
which  the  received  signal  level  due  to  auroral  scatter  was  esti¬ 
mated,  Cape  Romanzoff  on  the  coast  of  Alaska  was  selected. 
Environmental  conditions  at  the  primary  site  were  deemed  to  be 
unsuitable  for  the  equipment  and  a  secondary  location  in  Bethel 
in  southwest  Alaska  was  selected.  Figure  3-lb  shows  the  location 
of  the  Alaska  link. 
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Influenced  by  the  lessons  of  the  preliminary  Greenland 
experiment,  several  changes  to  the  basic  receive  equipment 
located  in  Bethel  were  made  (See  Figure  5-1).  A  slower  speed 
50  kHz/channel  A/D  converter  was  used  in  the  D-6000  data  acquisi¬ 
tion  system  due  to  the  reduced  sampling  requirements.  Data  was 
sampled  at  a  rate  of  5  samples  per  chip,  40  ksample s/sec . 
Instead  of  processing  data  immediately  after  acquisition  as  in 
Greenland,  raw  data  was  stored  on  flexible  disc  during  the  10 
minute  transmission  interval  and  processed  during  the  remaining 
50  minutes.  The  new  system  included  an  external  clock  with 
battery  backup  which  received  signals  from  WWV  to  allow  for  auto¬ 
matic  restart  after  power  failures.  The  experiment  was  designed 
to  run  for  six  months  with  only  bi-weekly  servicing  to  replace 
data  tapes.  The  equipment  in  Greenland  was  designed  for  manual 
restart  and  required  daily  service. 

The  low  bit  rate  and  potentially  high  Doppler  shifts  due  to 
scatter  from  the  rapidly  moving  aurora  limited  the  length  of  the 
probing  sequence  to  31  bits  (see  Equation  4.1  and  Table  5-1) 
although  a  longer  sequence  would  have  been  desirable  to  obtain 
improved  processing  gain. 

Post  acquisition  processing  on  the  data  which  was  acquired 
during  the  10  minute  probing  interval  was  performed  during  the 
remaining  50  minutes  of  each  hour.  During  the  first  four  minutes 
1024  sample  data  records  of  the  CW  signal  were  acquired  every 
five  seconds  and  1024  point  complex  FFTs  were  performed  on  every 
fifth  record.  The  magnitude  squared  of  the  FFT  record  was  com¬ 
puted  and  transferred  to  the  controller,  time  tagged,  and  stored 
on  tape.  Due  to  tape  storage  limitations,  phase  information  was 
not  transferred.  Cumulative  histograms  of  the  received  CW  signal 
level  were  performed  on  every  CW  record.  During  the  next  five 
minutes  of  the  10  minute  probing  period  the  PR  sequence  was 
repeated  continuously  by  the  transmitter.  A  data  window  corres¬ 
ponding  to  two  full  PR  sequence  periods  was  collected  every  six 

5-2 


.  — 4h  ft.  A  7  a  •  *.  *  _  .*  _  » 


Figure  5-1  Receive  Equipment  for  Alaska  Multipath  Measurement  Program 
Located  at  Bethel,  Alaska 


Table  5-1  Probing  Waveform  Characteristics 
for  Alaska  Experiment 


PR  SEQUENCE  CHARACTERISTICS 
Sequence  Length  (L) 

Maximum  Frequency  Offset 
Between  Rx  and  Tx 

Maximum  Doppler 

Shift  Expected 

Chip  Rate 

Chip  Period  (Tc) 

(Multipath  Resolution) 

Sequence  Repetition  Rate 
Maximum  Resolution 

Processing  Gain 
( Idealized ) 

Worst  Case  Expected. 

Processing  Gain  Degradation 

Number  of  Samples  Per- 
Sequence  Chip 

Sampling  Rate 

Snapshot  Rate 

F  FT  CHARACTERISTICS 

Repetition  Rate 

Number  of  Points 

F FT  Duration 

Frequency  Resolution 

F FT  Sampl  ing  Rate 

HISTOGRAM  CHARACTERISTICS 
Number  of  points 

Repetition  Rate  5 

Sampl  ing  Ra  te 
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31  bits 
±50  Hz 

±100  Hz 

8  kb/s 
12  5  us 

3.875  ms 

14.9  dB 

-2.6  dB 

5 

40  x  10^  samples/sec 
1  every  6  seconds 

1  every  25  seconds  (every 
5th  CW  record) 

1024  complex 

4  00  ms 

2.5  Hz 

390  us/sampie 

1024 

seconds — every  CW  record 
390  us/sample 


seconds  during  the  five  minute  interval.  Since  there  was  no 
timing  reference  to  align  the  data  acquisition  window  with  the 
beginning  of  the  PR  sequence,  two  full  PR  periods  were  acquired 
to  insure  that  one  full  contiguous  sequence  resided  within  the 
acquisition  window.  This  technique  greatly  simplified  the 
receive  equipment.  The  received  complex  data  window  was 
correlated  with  a  stored  version  of  the  PR  sequence.  The  magni¬ 
tude  squared  of  the  correlation  was  computed,  transferred  to  the 
controller,  time  tagged  and  stored  on  tape.  The  design  and  mix 
of  the  post-processing  routines  during  the  50  minute  idle  period 
were  driven  by  the  processing  capacity  of  the  Data  Precision  6000 
and  the  HP-85  controller  and  the  requirement  that  data  tapes  be 
changed  no  more  than  bi-weekly.  It  would  have  been  desirable  to 
acquire  more  FFT  records  and  to  save  both  magnitude  and  phase 
information;  however  the  limited  tape  capacity  precluded  this. 

During  the  construction  and  calibration  phase  of  the  exper¬ 
iment,  synthesizers .at  receiver  and  transmitter  were  measured  and 
determined  to  be  within  50  Hz  of  each  other.  Received  signal 
level  calibration  and  initial  noise  measurements  were  performed 
at  the  receive  site.  The  noise  measurements  showed  that  the 
receiver  at  Bethel  was  limited  by  man-made  noise  about  3  dB  above 
Galactic  noise.  Some  interference  due  to  images  from  the  several 
thousand  Watt  VOR  transmitter  at  the  nearby  airfield  were 
observed  about  200  KHz  from  the  operating  frequency.  They  were, 
for  the  most  part,  removed  by  the  30  KHz  IF  filter;  however,  this 
interference  represented  the  non-diurnal  noise  floor. 

5.2  RESULTS  FROM  ALASKA  EXPERIMENT 

Data  was  collected  over  a  seven  month  period  from  July  1985 
to  March  1986  with  several  interruptions  due  to  transmitter  or 
receiver  equipment  failures.  Due  to  the  remoteness  of  the 
receiver  location,  equipment  failures  often  required  four  weeks 
to  correct . 
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The  first  step  in  the  data  analysis  procedure  involved 
transferring  data  from  tape  cartridge  to  a  mainframe  computer 
and  removing  the  effect  of  tape  errors  which  often  caused  large 
spikes  in  the  data  records.  Next,  data  records  acquired  during 
periods  when  the  transmitter  was  not  operational  were  eliminated. 

Each  FFT  record  was  processed  to  determine  the  mean  and 
mean  squared  energy  of  the  data  record,  the  Doppler  shift  (mean), 
Doppler  spread  (second  central  moment),  and  the  skew  of  the 
Doppler  spectrum  (third  central  moment).  Cumulative  records  were 
formed  by  averaging  the  10  to  12  FFT  records  collected  during 
each  four  minute  acquisition  period.  The  various  moments  were 
computed  from  the  cumulative  FFT  record. 

Correlation  processing  was  more  complicated.  The  simplifi¬ 
cation  in  the  acquisition  process  which  eliminated  the  need  for 
aligning  the  acquisition  window  with  the  start  of  the  PR  sequence 
had  two  drawbacks  that  had  to  be  considered  in  the  processing. 
First  absolute  path  delay  could  not  be  computed  since  no  refer¬ 
ence  exists.  Second,  a  potential  ambiguity  in  the  magnitude  of 
the  delay  spread  had  to  be  removed  (see  Appendix  A).  Since  no 
time  reference  exists,  cumulative  averages  of  the  multipath  pro¬ 
file  could  not  be  produced  and  only  statistics  of  the  average  and 
variance  of  the  spread  measured  on  the  50  to  60  records  each 
five  minutes  were  computed. 

It  was  hoped  to  maintain  the  antennas  along  the  great 
circle  path  from  Anchorage  to  Bethel  to  observe  auroral  effects 
and  to  maximize  observation  of  sporadic-E  and  meteors.  Due  to 
low  solar  activity  during  the  period  of  the  experiments,  the 
aurora  remained  further  north  than  expected  and  several  changes 
to  the  experiment  plan  were  required.  We  divided  the  experiment 
into  two  phases . 

In  the  first  phase,  performed  from  mid-July  through  Octo¬ 
ber,  the  antennas  were  aligned  along  the  great  circle  path  from 
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Anchorage  to  Bethel  illuminating  the  midpoint  of  the  great  circle 
path.  Although  we  anticipated  the  observation  of  some  auroral 

scatter  activity,  the  primary  mechanism  that  it  was  hoped  to 
observe  was  meteors  and  daytime  sporadic-E.  In  the  second  phase 

of  the  experiment,  from  early  through  mid-July  and  from  late 

October  1985  through  March  1986,  the  antennas  were  rotated  off 
the  great  circle  path  to  observe  the  aurora  located  to  the 
North.  The  antenna  in  Anchorage  was  rotated  90  degrees  so  that 
it  pointed  in  the  direction  of  Fairbanks  and  the  antenna  at 
Bethel  was  rotated  45  degrees  north,  so  that  it  also  pointed 

towards  Fairbanks  (see  Figure  3-lb).  Rotating  the  transmitter 
antennas  required  collapsing  the  antenna  structure  to  allow 
manual  rotation  of  the  antennas.  Automatic  rotation  of  the 
antennas  to  track  the  aurora  was  therefore  not  possible. 

In  order  to  accurately  associate  multipath  conditions  with 
the  propagation  mechanisms  causing  them,  ionospheric  data  from 
the  University  of  Alaska  Ionosonde  at  College  (near  Fairbanks) 
and  various  geomagnetic  indicies  were  obtained. 

In  the  analysis  procedure  we  considered  the  characteristics 
of  propagation  mechanisms  which  were  similar  from  a  communication 
viewpoint  rather  from  an  ionospheric  physics  viewpoint.  For 
example,  we  grouped  together  all  forms  of  sporadic-E  although  on 
a  given  link  several  different  types  may  be  observed.  Also,  we 
separated  propagation  mechanisms  observed  on  and  off  the  great 
circle  path  from  transmitter  to  receiver. 

5.2.1  Sporadic-E  Propagation 

E-layer  propagation  including  enhanced  ionospheric  scatter 
and  sporadic-E  propagation  were  observed — though  not  as  fre¬ 
quently  as  on  some  of  the  longer  links  in  Alaska.  For  the  pur¬ 
pose  of  communication  analysis  we  group  together  all  forms  of 
sporadic-E  and  enhanced  E-layer  ionospheric  scatter  under  the 
general  heading  sporadic-E. 
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The  sporadic-E  channel  was  characterized  by  signal 
strengths  in  excess  of  8  dB  above  the  noise  floor  in  the  30  KHz 
bandwidth,  slow  fading  (2o  fading  bandwidth  usually  less  than 
20  Hz)  and  little,  if  any,  multipath  spreading  within  the  125  us 
minimum  resolution  of  the  experiment.  The  primary  multipath 
mechanism  observed  was  the  occasional  occurrence  of  meteors  in 
addition  to  the  sporadic-E  signal. 

Figure  5-2  shows  a  series  of  multipath  snapshots  of  a 
typical  sporadic-E  event  acquired  every  six  seconds  for  five 
minutes  on  1  July  1985  at  0300  UT.  The  sporadic-E  layer  initial 
frequency,  fgEs,  measured  at  College  Alaska  was  5.14  MHz.  The 
frames  (snapshots)  have  been  aligned  so  that  all  the  peaks  occur 
at  the  same  relative  delay.  In  snapshots  21,  33,  38,  44,  and  49 
the  secondary  discrete  components  are  due  most  likely  to  meteor 
trails  to  the  side  of  the  mid-point  of  the  great  circle  plane. 
The  delays  correspond  to  path  differences  of  30,  112,  75,  37,  and 
75  Km  respectively.  Figure  5-3  shows  the  FFT  of  the  CW  received 
probing  signal  averaged  over  the  4  minutes  preceding  the  multi- 
path  snapshots.  The  measured  2a  Doppler  spread  is  11.8  Hz. 

Since  the  time  between  multipath  snapshots  (6  seconds)  is 
much  greater  than  the  inverse  of  the  fading  bandwidth  in  Fig¬ 
ure  5-3,  it  is  expected  that  the  squares  of  the  correlation  peaks 
minus  the  noise  floor  (signal  intensity  I)  represent  the  energy 
in  the  signal  envelope  and  should  be  uncorrelated  from  sample  to 
sample.  If  the  fading  of  the  correlation  peaks  were  due  to 
multiple,  closely  spaced  sporadic-E  echoes,  then  one  would  expect 
the  received  signal  to  have  complex  Gaussian  statistics  and  an 
exponential  distribution  of  signal  energies.  Figure  5-4  shows 
the  histogram  of  the  square  of  the  peaks  minus  the  noise  floor. 
The  broken  line  shows  a  MMSE  exponential  curve  fit  to  the  data. 
We  observe  a  relatively  close  fit. 
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Another  statistical  parameter  that  was  computed  was  the 
Scintillation  Index,  Sj,  defined  as 


where  I  is  the  signal  intensity  (square  of  the  envelope  with  the 
noise  floor  removed.)  For  the  data  of  Figure  5-2,  with  histogram 
shown  in  Figure  5-4,  the  Scintillation  Index  computed  from 
Equation  (5.2)  was  0.89.  The  average  Scintillation  Index  for  all 
sporadic-E  events  was  1.21  with  a  standard  deviation  of  0.4. 
Scintillation  Index  equal  to  zero  corresponds  to  a  constant  non¬ 
fading  signal,  one  to  a  fading  signal  with  complex  Gaussian 
statistics,  and  indices  larger  than  one  correspond  to  data  sets 
in  which  the  number  of  multipath  components  is  small  with  one 
component  larger  than  the  others. 

Figure  5-5  shows  a  histogram  of  the  2a  Doppler  spreads  of 
Sporadic-E  propagation  data  obtained  with  the  antennas  on  the 
great  circle  path.  The  average  Doppler  spread  over  ail 
Sporadic-E  events  was  10.958  Hz  with  a  3.25  Hz  standard  devia¬ 
tion.  The  10,50,  and  90  percentile  Doppler  spreads  were  7  Hz, 
10  Hz,  and  17  Hz,  respectively. 

Figure  5-6  shows  the  percentage  of  time  sporadic-E  was 
observed  vs  time  of  day.  Sporadic-E  propagation,  as  shown  in 
Figure  5-6,  was  observed  most  frequently  around  local  noon 
(2100  CJT)  with  the  antennas  pointed  along  the  great  circle 
path.  Sporadic-E  propagation  was  not  observed  with  the  antennas 
pointing  towards  Fairbanks  since  the  region  in  the  vicinity  of 
the  mid-point  of  the  great  circle  path  was  not  illuminated  by  the 
main  lobe  of  the  antenna  patterns.  Table  5-2  summarizes  the 
characteristics  of  the  sporadic-E  channel. 

We  now  examine  briefly  the  relationship  between  observed 
sporadic-E  conditions  and  ionosonde  parameters.  The  primary 


Table  5-2  Characteristics  of  Sporadic-E  Channel 


Number  of  Observation  Hours 

Number  of  Observation  Days 

Number  of  Sporadic-E  Hours 

Number  of  Sporadic-E  Days 

Average  2o  Doppler  Spread 

Standard  Deviation 

10  Percentile  Doppler  Spread 

50  Percentile  Doppler  Spread 

90  Percentile  Doppler  Spread 

Multipath  on  Sporadic-E  Channel 

Average  Scintillation  Index 

Standard  Deviation  of  Scintillation  Index 

Average  SNR  (Sporadic-E) 

30  KHz  BW 


738 

45 

106  (14.3%) 
20  (44.44%) 
10.95  Hz 
3.25  Hz 
7  Hz 
10  Hz 
17  Hz 
Not  Measurable 
1.21 
0.4 
8.9  dB 
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measure  available  was  the  E-layer  and  sporadic-E  critical  fre¬ 
quencies  f q E  and  fgEs  determined  by  verticle  swept  sounding  at 
College,  Alaska.  Figure  5-7  shows  the  correspondence  between  the 
maximum  of  f0E  and  fgEs  and  the  total  received  signal  power.  The 
broken  line  represents  maximum  ( f0E  and  fgEs)  at  College,  Alaska 
at  zero  minutes  past  each  hour.  The  solid  line  represents  the 
total  received  signal  power  minus  the  noise  floor  in  uncalibrated 
power  units  determined  by  integrating  the  energy  in  the  composite 
FFT .  The  days  shown  are  30  June  through  2  July  1985  with  the 
antennas  on  the  great  circle  path. 

Note  that  with  the  exception  of  one  measurement  in  which  a 
very  high  fgEs  was  recorded  at  College  and  little  sporadic-E 
propagation  was  observed,  the  occurrence  of  sporadic-E  propaga¬ 
tion  on  the  Anchorage-Bethel  link  closely  tracks  the  fgE  and  foEs 
at  College. 
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From  a  commun  ications  point  of  view,  sporadic-E  propagation 
results  in  a  continuous  channel  with  high  signal  strength,  slow 
fading  and  low  multipath  spread,  providing  enhanced  throughput 
capacity  for  meteor  burst  systems.  Ironically,  the  occasional 
occurrence  of  meteors  represents  the  primary  source  of  multipath 
interference  when  the  sporadic-E  channel  is  present.  Sporadic-E 
propagation  can  support  coherent  modulation  techniques  with 
coherence  bandwidths  often  on  the  order  of  Megahertz.  Link 
protocols  such  as  ARQ  (automatic  request)  or  Acknowledge/Non- 
acknowledge  are  adequate  for  communications  although  they  are  not 
optimal  for  the  slow  fading  channel.  Diversity  reception  may  be 
necessary  to  fully  exploit  the  bandwidth  available.  Advanced 
techniques  which  can  adapt  the  data  rates  to  the  slowly  changing 
capacity  of  the  channel  could  be  used  to  fully  exploit  the  chan¬ 
nel  resource  [Weitzen,  1983]. 
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At  the  network  level,  protocols  must  be  capable  of  full  and 
normal  operation  in  the  high  contention  environment  created  by 
sporadic-E  propagation  in  which  continuous  connectivity  to  many 
nodes  in  a  network  may  exist  for  extended  periods  of  time.  Pure 
meteor  scatter  network  protocols  operate  on  the  principle  that 
the  duty  cycle  on  any  given  link  is  so  low  that  collisions  be¬ 
tween  network  transmissions  are  negligible.  A  pure  meteor  scat¬ 
ter  link  could  become  collision  bound  in  a  frequent  sporadic-E 
env  ironment . 

The  primary  drawback  with  the  sporadic-E  channel  is  that  it 
is  not  as  predictable  or  reliable  as  meteors.  It  is  much  more 
subject  to  ionospheric  conditions  than  meteors.  Therefore  from  a 
reliablity  stand,  it  should  be  treated  as  a  "bonus"  channel 
rather  than  as  the  primary  channel  when  it  is  used  as  a  surviv- 
able  backup  communication  technique.  Because  the  sporadic-E 
channel  has  a  much  larger  scattering  area  than  meteors,  it  has 
much  less  resistance  to  intentional  or  unintentional  interference 
than  do  meteors . 

5.2.2  Multipath  Due  To  Multiple  Meteors 

While  not  specifically  designed  for  the  observation  of 
meteor  trails  (due  to  the  fact  that  the  receiver  processor  trig¬ 
gers  at  fixed  time  intervals  rather  than  on  intermittent  signals 
such  as  meteors)  ,  a  number  of  meteors  were  observed  on  the 
link.  Meteors  were  observed  with  the  antennas  on  the  great 
circle  path  and  pointing  north,  although  as  expected,  many  more 
meteors  were  observed  with  the  antennas  on  the  great  circle  path. 

On  occasion,  multipath  effects  due  to  the  simultaneous 
occurrence  of  more  than  one  meteor  trail  were  observed.  Using 
data  sets  with  the  antennas  on  the  great  circle  path  with 
sporadic-E  effects  removed,  a  six  percent  duty  cycle  for  meteor 
trails  was  observed  with  the  10,000  watt  transmitter.  Of  the  six 


5-18 


percent  of  all  records  which  had  meteor  trails,  1.5  percent  (of 
the  six  percent)  showed  secondary  components  due  to  multiple 
meteor  trails.  This  compares  favorably  to  the  one  to  three 
percent  multipath  predicted  by  Sugar  [I960]. 

5.2.3  Auroral  Scatter 

In  regions  in  the  vicinity  of  the  aurora,  scatter  from  the 
aurora  represents  a  serious  potential  impairment  to  meteor 
communication  at  high  latitudes.  Scatter  from  the  aurora  can 
transform  the  relatively  benign,  though  intermittent,  meteor 
channel  with  little  multipath  interference  into  a  continuous 
fading  channel  with  2  o  delay  spreads  in  excess  of  1000  us  and  2a 
Doppler  bandwidths  in  excess  of  200  Hz. 

Auroral  measurements  were  conducted  with  the  antennas 
pointing  north  from  early  to  mid-July  1985  and  from  late  October 
through  March  1986. 

Auroral  effects  were  observed  most  frequently  but  not 
exclusively  during  the  late  night  and  early  morning  hours, 
11  P.M.  to  6  A.M.  local  time  (0800  to  1500  UT )  and  the  observa¬ 
tion  of  auroral  effects  could  be  correlated  to  increases  in  the 
standard  geomagnetic  planetary  index  Kp .  Figure  5-8  shows  a 
sequence  of  multipath  snapshots  collected  every  six  seconds 
beginning  at  1205  UT  on  1  November  1985.  Figure  5-9  shows  the 
composite  Doppler  profile  formed  by  averaging  10  individual  FFT 
profiles  acquired  every  25  seconds.  Figures  5-10  and  5-11  show 
multipath  sequences  and  the  Doppler  profile  respectively  for  a 
period  two  hours  later.  Figures  5-12  and  5-13  show  the  multipath 
and  Doppler  profiles  for  the  spectacular  event  on  the  evening  of 
30  November  1985  .  Figure  5-14  and  5-15  show  the  multipath 
Doppler  profiles  for  a  typical  auroral  event  on  July  4,  1985. 

In  Figures  5-8,  5-10,  5-12,  and  5-14,  the  noise  floor  has 

been  removed  for  clarity  and  the  multipath  profiles  have  been 
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aligned  manually  so  that  components  appear  in  the  same  time  bins 
(recall  that  there  is  no  absolute  time  reference).  Observe  in 
the  multipath  sequences  how  components  appear,  fade  and  then 
disappear  independently,  implying  an  uncorrelated  scatter  fading 
channel . 

Figures  5-8  through  5-15  show  the  characteristics  of  the 

purely  auroral  channel  observed  with  the  antennas  off  the  great 

circle  path.  Because  the  antennas  were  pointed  far  off  the  great 
circle  path,  few  meteors  were  observed  and  no  sporadic-E  was 

observed.  The  2 o  Doppler  and  multipath  spreads  in  Figures  5-8 
and  5-9  are  210  Hz  and  589  ms,  respectively.  In  Figures  5-10  and 
5-11  the  Doppler  and  multipath  spreads  are  190  Hz  and  624  us, 

respectively,  and  in  Figures  5-14  and  5-15  the  Doppler  and 
multipath  spreads  are  110  Hz  and  825  ms,  respectively.  For  the 
spectacular  display  of  30  November,  1985,  the  Doppler  and  multi- 
path  spreads  were  246  Hz  and  1005  ms,  respectively.  Many  auroral 
events  were  faint  so  that  they  were  apparent  from  the  FFT  frames 
(400  ms  records,  effective  processing  gain  *  30  dB)  and  not 
visible  or  faintly  visible  on  the  multipath  records  with  less 
than  15  dB  processing  gain. 

We  now  determine  several  comprehensive  statistics  for  the 
pure  auroral  scatter  channel.  Figure  5-16  shows  a  histogram  of 
2a  Doppler  spreads  for  all  events  classified  as  auroral  scatter 
with  the  antennas  off  the  great  circle  path.  The  average  2a 
Doppler  spread  for  all  events  was  100.3  Hz  with  a  67  Hz  standard 
deviation.  The  10,50,  and  90  percentile  Doppler  spreads  were 
47  Hz,  105  Hz,  and  200  Hz,  respectively.  Figure  5-17  shows  the 
histogram  of  2o  multipath  spreads  for  each  event  identified  as 
auroral  scatter.  The  average  multipath  spread  was  804  us  with  a 
174  us  standard  deviation.  The  10,  50,  and  90  percentile  multi- 
path  spreads  were  420  us,  800  us,  and  1000  us,  respectively. 
Table  5-3  summarizes  the  data  in  terms  of  basic  communication 
parameters . 
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Table  5-3  Auroral  Scatter  Results 


Number  of  Observation  Hours 

Number  of  Hours  Auroral 
Scatter  Observed 

Average  2c  Doppler  Spread 

STD  Deviation  of  Doppler  Spread 

10  Percentile  Doppler  Spread 

50  Percentile  Doppler  Spread 

90  Percentile  Doppler  Spread 

Average  2o  Multipath  Spread 

STD  Deviation  of  Multipath  Spread 

10  Percentile  Multipath  Spread 

50  Percentile  Multipath  Spread 

90  Percentile  Multipath  Spread 

Average  SNR  During  (1  kHz  BW) 
Auroral  Events 


74  5 

34  (4.5%) 

100.3  Hz 
67  Hz 
47  Hz 
105  Hz 
200  Hz 
80  6  us 
174.23  us 
42  0  us 
80  0  us 
100  0  us 
13  dB 


As  observed  in  Figures  5-8,  5-10,  5-12,  and  5-14,  received 
auroral  scatter  signals  generally  consist  of  the  sum  of  a  number 
of  independently  fading  components.  The  envelope  squared  of  the 
signal  should  therefore  be  exponentially  distributed.  Fig¬ 
ure  5-18  shows  a  histogram  of  the  received  envelope  squared  for 
the  time  period  shown  in  Figures  5-14  and  5-15.  The  broken  line 
represents  a  MMSE  exponential  fit  to  the  data.  We  see  that  there 
appears  to  be  a  good  fit  to  the  data.  Due  to  a  calibration 
error,  absolute  signal  levels  are  not  available  in  Figure  5-18. 

The  night  of  September  1  9,  1985  provided  us  with  an  oppor¬ 
tunity  to  observe  aurora  with  the  antennas  on  the  great  circle 
path  to  show  what  a  meteor  channel  would  look  like  with  high 
level  background  aurora.  Figure  5-19  shows  the  sequence  of 
multipath  snapshots  for  the  September  19  event.  Large  meteors 
above  the  background  aurora  are  observed  in  frames  47,  53,  and 
54.  Figure  5-20  shows  the  Doppler  profile  in  which  we  can  see 
two  distinct  components;  the  zero  shift  component  due  to  the 
meteor  trails,  and  the  more  diffuse  component  centered  at  about 
80  Hz  represents  the  auroral  scatter. 

The  observation  of  Auroral  Effects  on  the  Anchorage-8e thel 
link  shows  a  strong  correlation  to  increases  in  the  geomagnetic 
index  Kp.  Figure  5-21  shows  the  received  signal  energy  due  to 
scatter  from  the  aurora  (solid  line)  in  uncalibrated  signal  units 
and  the  three  hour  geomagnetic  index  Kp  for  4  to  9  July  1985 
(broken  line)  highlighted  in  Figures  5-14  and  5-15.  -Note  that 
there  appears  to  be  several  hours  delay  between  the  peak  in  Kp 
and  the  observation  of  auroral  scatter  effects. 

Figure  5-22  shows  the  received  signal  energy  due  to  auroral 
scatter  and  the  Kp  magnetic  index  for  the  period  1  November 
through  4  November  1985  ,  highlighted  in  Figures  5-8  through 
5-11.  Figure  5-23  shows  the  received  signal  energy  due  to 
auroral  scatter  and  the  Kp  magnetic  index  (broken  line)  for  the 
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Figure  5-20  Composite  Doppler  Profile  for  Period  Shown  in  Figure  5-19. 

Component  Centered  on  Zero  Shift  Corresponds  to  Meteor  Propagation 
and  the  Shifted  Component  is  Due  to  the  Auroral  Scatter. 
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Figure  5-22  Solid  Line  Represents  Received  Signal  Energy  in  Uncalibrated  Signal 
Units  for  1  to  3  November  1985.  Broken  Line  Represents  3  hour  Kp 
Geomagnetic  Activity  Index. 


period  29  November  through  5  December  1985,  highlighted  in  Figure 
5-12  and  5-14.  Figure  5-24  shows  the  received  signal  energy  to 
auroral  scatter  (solid  line)  and  Kp  magnetic  index  (broken  line) 
for  the  period  17  September  through  21  September  1985, 

highlighted  in  Figures  5-19  and  5-20. 

The  observation  of  auroral  effects  could  also  be  correlated 
to  the  absence  of  fgE  and  fgEs  traces  on  the  ionosonde  data. 
This  is  attributable  to  the  fact  the  increased  ionization  levels 
absorb  the  low  frequency  ionosonde  probing  signal. 

Figures  5-22,  5-23,  and  5-24  show  the  Kp  index  (broken 

line)  and  received  signal  level  (solid  line)  for  the  1  to 

3  November,  29  November  to  3  December  1985,  and  September  17  to 

20  periods  highlighted  in  Figures  5-8  through  5-13,  5-19,  and 

5-20.  Finally,  Figure  5-25  shows  the  percentage  of  days  on  which 
auroral  scatter  effects  were  observed  versus  time  of  day  for  the 
entire  period  of  the  experiment  with  the  antennas  off  the  great 
circle  path.  Note  that  the  hour  between  19  and  00  hours  corres¬ 
ponding  to  local  noon  shows  the  lowest  activity. 


5.2.3.  1 


Malaga  [1986  ]  has  developed  in  detail  the  theory  of  VHF 
scatter  from  the  aurora.  The  average  power  received  due  to  scat¬ 
ter  from  the  aurora  is  determined  by  Malaga  as 
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where 


f()E  is  the  E-layer  critical  frequency 
PT 
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is  the  transmitter  power 


age  of  Days  on  which  Auroral  Scatter  Effects  Were  Observed 
Time  of  Day  ( UT ) .  Local  Daylight  Time  is  UT-9  Hours. 


Gt  is  the  peak  transmitter  antenna  gain 
Gr  is  the  peak  receiver  antenna  gain 
c  is  the  speed  of  light 

f  is  the  carrier  frequency 

%  is  the  two-way  absorption  loss  incurred  as  the 
incident  and  scattered  waves  propagate  through 
the  ionosphere. 

CE(m)  =  (m-2)L  (— is  the  structure  parameter  of 

the  irregularities 
N  is  the  mean  electron  density 

AN  is  the  perturbation  in  electron  density  due  to  field 
aligned  irregularities. 

L  is  the  outer  scale  of  the  irregularities  along  the 
magnetic  field  lines. 

Lf  is  the  outer  scale  of  the  irregularities  in  the  plane 
transverse  to  the  magnetic  field. 

m  is  the  wavenumber  spectrum  slope  which  determines  the 
frequency  dependence  of  the  scattered  field. 
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integrates  the  contribution  of  the  scatterers  over  the 
three  dimensional  antenna  pattern. 

g^ ( h , 6 , $ )  is  the  transmitter  antenna  pattern 

gr(h,0,$)  is  the  receiver  antenna  pattern 

R-jvR^  are  the  distances  to  the  scatterer  in  the  common 
volume  from  transmitter  and  receiver,  respectively. 

uf  (h,0,$)  is  a  unit  vector  parallel  to  the  direction  of 
the  magnetic  field 

ui  (h,0,$)  is  the  unit  vector  directed  from  transmitter 
to  the  scatterer. 

Ug  (h,0,<j>)  is  the  unit  vector  directed  from  the  scatterer 
to  the  receiver. 

VL ( h  ,  0 ,  <J> )  *  (uQ-ui)  •  uf 

V  (  h  ,  0 ,  <p  )  =  lug-uj2 

VT(h,0,$)  =  v2~V2 

X  is  the  angle  between  incident  E  vector  and  direction 
of  scattering. 

Kl  =  ^  ^ 

Equation  (5.2)  assumes  that  the  critical  frequency  of  the 
scattering  layer,  fgE,  and  the  normal  electron  density 
fluctuations  are  constant  within  the  region.  The  triple  integral 
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I  (h0,b)  depends  only  on  the  geometry  and  polarization  of  the 
transmit  and  receive  antennas  and  can  be  evaluated  given  the 
geographic  location  of  the  terminals  and  the  scattering  region. 
Prediction  of  the  scattered  power  requires  modeling  or  measuring 
the  critical  frequency  fQE»  the  absorption  loss  K£,  and  the 
structure  parameter  CE.  Due  to  the  geometry  of  the  scattering 
(see  Malaga  [1986]')  ,  signals  observed  at  the  earth  are  generally 
due  to  scatterers  in  the  E  region  of  the  ionosphere. 

The  rms  delay  spread  can  be  calculated  from  Equation  (5.2) 
by  computing  various  moments  of  the  delay  scaled  by  the  power 
from  all  scatterers  in  the  common  volume  as 


P  is  the  received  power  from  each  scatter  in  the  common 
_  volume  as  a  function  of  location 

kT 

— - —  is  the  path  delay  from  transmitter  to  scatterer  to 
receiver  at  each  point  in  the  common  volume. 


The  2a  delay  spread  is  given  by 


Define  the  total  path  loss  Lp  as 


L  =  •  r-V  =  — 2 - - -  •  <  5  •  5  ) 

*T  T  R  Cg  exp  (-0.23kg) 

Lr  is  the  relative  path  loss  which  accounts  for  all  propagation 

effects  except  absorption  loss,  Kg,  and  terms  incorporated  in  the 

2 

structure  parameter  C£. 

Figure  5-26  shows  contours  of  relative  path  loss  versus 
receiver  location  in  Alaska  calculated  from  Equation  (5.2)  and 
(5.3)  assuming  a  42.4  MHz  transmitter  in  Anchorage,  Kp  =  2,  the 
auroral  region  bounds  due  to  Gassmann  [1973],  and  10  km  longi¬ 
tudinal  outer  scale  of  irregularities.  Figure  5-27  shows 
contours  of  delay  spread  versus  receiver  location  in  Alaska  with 
a  transmitter  located  in  Anchorage.  The  locations  of  Anchorage 
and  Bethel  are  shown  on  the  figures. 

Note  that  computation  of  delay  spread  in  Figure  5-27  does 
not  require  knowledge  of  the  physical  parameters  incorporated  in 
the  structure  parameter  or  the  absorption  loss.  For  Anchorage- 
Bethel  we  compare  the  1000  ms  average  rms  delay  spread  prediction 
to  the  806  us  delay  spread  observed  on  the  link.  We  see  from 
Figure  5-26  that  with  a  transmitter  in  Anchorage,  Bethel  would 
appear  to  be  a  good  location  for  the  equipment. 

Accurate  comparison  of  the  absolute  received  signal  predic¬ 
tions  to  data  would  require  knowledge  of  the  absorption  KE ,  the 

2 

E-layer  critical  frequency  fgE  and  the  structure  constant  Cp,  or 
given  the  data,  f0E  and  KE,  we  could  calculate  C£. 

In  order  to  acquire  this  data,  receive  and  transmitter 
equipment  would  require  frequent  and  accurate  calibration.  In 
addition,  accurate  antenna  pattern  measurements  would  be  re¬ 
quired.  Simple  f g E  measurements  using  the  College  vertical  swept 
frequency  ionosonde  equipment  would  not  have  been  adequate  due  to 
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Figure  5-27  Contours  of  constant  delay  spread  in  microseconds 


the  fact  that  the  low  frequency  probing  waveform  was  often 
absorbed.  An  oblique  angle  prober  would  have  been  required.  Due 
to  the  costs  required  to  accurately  determine  these  parameters, 
it  was  decided  to  concentrate  on  characterizing  the  multipath  and 
Doppler  spreads  which  could  be  determined  without  the  expensive 
calibration  and  using  existing  ionospheric  probers. 

5.2.3. 2  The  Auroral  Scatter  Communication  Channel 

Auroral  scatter  at  low  levels  represents  a  potential  back¬ 
ground  noise  source.  At  high  signal  levels,  measurements 
described  in  this  section  show  it  to  represent  a  highly  disturbed 
channel  with  multipath  spreads  often  in  excess  of  1000  us  and 
Doppler  spreads  of  200  Hz  or  more  as  shown  in  Figures  5-19  and 
5-20.  When,  and  if,  the  auroral  channel  is  dominant,  convention¬ 
al  meteor  burst  communication  at  4-8  Kbps  using  coherent  modula¬ 
tion  is  virtually  impossible.  Even  non-coherent  communication  at 
these  rates  would  be  intersymbol  interference  limited. 

At  a  minimum,  protocols  for  a  meteor  burst  network  opera¬ 
ting  in  an  auroral  environment  must  be  designed  so  that  one  link 
operating  on  a  continuous  though  unusable  channel  does  not  jam  or 
lockout  the  rest  of  a  network  which  may  not  be  experiencing 
auroral  scatter  effects.  A  lower  rate  (100  to  300  bps)  non¬ 
coherent  operating  mode  might  be  included  in  a  high  latitude 
system  just  in  case  the  link  becomes  auroral  scatter  limited. 

Disruption  of  meteor  communication  due  to  scatter  from  the 
aurora  is  a  localized  condition  which  depends  on  the  antenna 
patterns  of  a  link  and  on  the  location  and  orientation  of  the 
aurora  with  respect  to  the  link  and  the  antenna  patterns.  Not 
all  links  in  the  vicinity  of  the  aurora  will  be  effected.  The 
Bethel-Anchorage  link  was  selected  for  the  experiment  because  it 
was  predicted  to  be  one  of  the  worst  (or  best  depending  on  how  we 
look  at  it)  links  for  observing  auroral  maultipath  effects  given 


a  transmitter  in  Anchorage.  Some  links  may  experience  some 
slight  effects  while  other  links  may  satisfy  the  requirements  for 
scatter  from  the  aurora  and  may  experience  the  disturbed  scatter 
conditions.  Malaga  [1986]  describes  the  exact  conditions  re¬ 
quired  for  auroral  scatter  and  presents  some  predictions  for 
scatter  conditions  in  Alaska. 

To  actually  communicate  over  the  auroral  channel,  advanced 
adaptive  equalization  techniques  are  required  such  as  those 
applied  to  troposcatter ,  or  underwater  communication  systems. 
The  alternative  is  to  lower  the  data  rate  well  below  the 
multipath  spread  and  operate  using  non-coherent  modulation  and 
cod ing . 
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SECTION  6 

SUMMARY:  COMMUNICATING  VIA  METEOR  BURST  AT  HIGH  LATITUDES 


A  number  of  propagation  mechanisms  in  addition  to  meteors 
are  observed  on  the  high  latitude  meteor  channel.  In  this  report 
we  have  characterized  the  multipath  and  fading  profiles  of  the 
dominant  mechanisms:  meteors,  sporadic-E,  and  auroral  scatter. 

We  have  observed  that  some  of  the  mechanisms  can,  if  link  proto¬ 
cols  are  designed  properly,  be  beneficial. 

High  latitude  sporadic-E  and  enhanced  E-layer  scatter  can 
create  a  continuous,  slow  fading,  low  multipath  communication 
channel  similar  to  a  very  benign  single  layer  HF  channel.  In 
order  to  exploit  this  channel,  network  protocols  must  be  care¬ 
fully  designed  so  that  the  network  can  operate  in  a  high  conten¬ 
tion  environment  when  all  or  many  links  have  continuous  connec¬ 
tivity  to  the  master  station.  Variable  rate  modems  capable  of 
tracking  the  slow  fading  channel  would'  allow  optimal  use  of  the 
capacity  of  the  channel .  Sporadic-E  channels  have  been  observed 
on  many  high  latitude  links  and  often  confused  for  meteor 
channels.  The  sporadic-E  channels  do  not  have  the  inherent 
AJ/LPI  of  meteor  channels  since  they  have  a  much  larger  scatter¬ 
ing  area  than  meteors.  They  are  also  not  as  reliable  as  is  the 
meteor  channel  and  are  much  more  subject  to  ionospheric  disturb¬ 
ances  than  is  the  meteor  channel.  It  therefore  should  be  treated 
as  a  "bonus"  channel  rather  than  the  backbone  for  survivable/ 
reliable  austere  backup  communication.  If  the  meteor  channel 
alone  cannot  support  the  throughput  required  for  the  mission, 
depending  on  the  sporadic-E  channel  for  the  throughput  require¬ 
ments  is  potentially  dangerous. 

Scatter  from  the  aurora  can  produce  a  low  level  background 
interference  source  on  some  links  and  depending  on  solar  activity 
and  the  location  of  the  aurora  relative  to  a  link  could  disrupt 
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meteor  communication  on  some  links.  When  dominant,  the  auroral 
channel  is  a  highly  disturbed  channel  with  multipath  spreads 
potentially  in  excess  of  1000  us  and  Doppler  spreads  in  excess  of 
200  Hz.  Communication  using  conventional  8  kbps  coherent  BPSK 
modulation  on  such  a  channel  would  be  virtually  impossible.  Low 
rate  ( i  .e .  )  100-300  bps  non-coherent  modulation  with  frequency¬ 
hopping  and  coding  would  be  better  suited  to  the  channel  unless 
advanced  equalization  techniques  are  incorporated.  With  the 
lower  data  rate,  a  continuous  channel  with  a  throughput  on  the 
same  order  as  the  intermittent  higher  burst  rate  meteor  channel 
could  be  achieved.  Communication  protocols  must  be  very 

carefully  designed  so  that  a  link  with  a  continuous,  though 
unusable,  channel  to  the  master  station  does  not  disrupt  or 
lockout  the  entire  network. 

Taking  into  account  the  undesirable  occasional  effect  of 
auroral  scatter,  meteor  burst  still  has  many  advantages  over  HF 
which  can  be  even  more  effected  by  these  same  phenomena.  Auroral 
scatter  degradation  at  meteor  burst  frequencies  may  be  absorption 
blackout  conditions  at  HF . 

The  experimental  program  as  a  whole  was  a  success.  The 
major  goals  of  characterizing  the  multipath  and  Doppler  profiles 
of  the  meteor,  sporadic-E  and  auroral  Channels  with  a  short  dur¬ 
ation,  relatively  inexpensive,  experiment  were  satisfied. 

Interference  with  sensitive  air  traffic  control  equipment 
at  Sondestrom  Air  Base,  Greenland  by  the  wideband  probing  signal 
forced  early  termination  of  the  wideband  preliminary  experiment 
in  Greenland.  The  results  from  that  experiment  were  somewhat 
inconclusive  due  to  the  limited  data  set;  however,  no  multipath, 
other  than  that  caused  by  meteor  trails  was  observed  on 
sporadic-E  channels  within  the  8  us  resolution  of  the  experi¬ 
ment.  The  lessons  learned  during  the  experiment  were  applied  to 
the  main  experiment  in  Alaska  and  contributed  to  its  success. 


The  three  primary  propagation  mechanisms  observed  at  high 
latitudes,  meteors,  sporadic-E  ,  and  auroral  scatter  were  ob¬ 
served  and  characterized  by  the  Alaska  experiment. 

6.1  RECOMMENDATIONS  FOR  FUTURE  STUDY 

The  work  described  in  this  and  the  companion  report 
[Malaga,  1986]  represent  an  important  first  step  in  scientif¬ 
ically  evaluating  the  reliability  of  meteor  burst  at  high  lati¬ 
tudes.  The  research  results  must  now  be  applied  to  practical 
communications  problems.  A  new  meteor  burst  communication  net¬ 
work  will  become  operational  in  1986  to  link  remote  long  range 
radar  (LRR)  sites  throughout  Alaska  to  a  central  facility  near 
Anchorage.  The  network  will  provide  austere  backup  communication 
to  the  normal  satellite  communication  systems.  Figure  6-1  shows 
the  location  of  the  meteor  burst  sites  in  the  network.  In 
Figure  6-2  we  superimpose  a  map  showing  the  location  of  the  radar 
sites  upon  the  contour  maps  of  relative  path  loss  computed  by 
Malaga.  The  figure  shows  which  sites  might  be  prime  candidates 
for  auroral  degradation. 

Based  on  these  predictions  from  the  Malaga  report,  the  fol¬ 
lowing  links  in  the  new  Alaska  Air  Command  system  would  be  prime 
candidates  for  performance  degradation  due  to  auroral  scatter. 

Anchorage-Cape  Romanzoff 
Anchorage-Tatal  ina 
Anchorage-Cape  Newenham 
Anchorage-Sparrevohn 

In  addition,  an  Anchorage-Fa irbanks  link  would  be  especial¬ 
ly  susceptible-  During  disturbed  conditions,  more  southerly  sta¬ 
tions  might  be  effected. 


6-3 


'Cape  Llsbume 


LKotzebue 


'Tin  City. 


•  Fort  Yukon 
*  Indian  Mountain  ^ 

•  Galena  #  Murphy  Dome 


1  Cape  Romanzof 


•  TataJlna 


•  Sparrevohn 


Cape  Newenham 


King  Salmon 


Anchorage 
(master  station) 


3  o'  rs 

Cold  Bay 


Figure  6-1  Location  of  Meteor  Burst  Stations  in  Alaskan 
Air  Command  Network 


A  useful  and  valuable  follow-on  to  the  experiment  described 
in  this  report  would  be  to  collect  data  from  the  Alaskan  Air  Com¬ 
mand  on  the  performance  of  the  various  links  in  the  network  to 
assess  which  links  are  dominated  by  sporadic-E  channels,  which 
links  are  primarily  meteor  channels  and  which  links,  if  any,  are 
limited  by  auroral  scatter  effects. 

On  the  links  showing  auroral  scatter  performance  degrada¬ 
tion,  instrumentation  such  as  used  in  this  experiment  should  be 
set  up  to  determine  the  multipath  spreads  on  the  links  and  to 
correlate  measured  multipath  spreads  with  communication  perfor¬ 
mance  degradation.  If  auroral  scatter  proves  to  be  a  problem  on 
the  network,  techniques  to  mitigate  the  effects  including  adap¬ 
tive  antenna  steering,  increased  operating  frequencies,  reduced 
data  rates,  or  adaptive  message  routing  should  be  explored. 
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APPENDIX  A 

REMOVING  THE  AMBIGUITY  IN  THE  MULTIPATH  SPREAD 


To  show  how  an  ambiguity  in  the  measured  delay  spread  can 
be  created  by  the  data  processing  algorithms  and  then  removed, 
consider  a  multipath  profile  that  consists  of  two  discrete 
components  with  relative  delay  0.75  P,  where  P  is  the  period  of 
the  PR  probing  sequence.  The  location  of  the  multipath  profile 
relative  to  the  data  acquisition  window  is  random  from  window  to 
window.  Consider  first  the  case  in  which  the  first  component  is 
located  at  the  beginning  of  the  data  acquisition  window  and  the 
second  component  is  delayed  by  0.75  P.  In  this  case,  there  is 
no  ambiguity.  Next  consider  the  case  in  which  the  window  is 
aligned  so  that  the  first  component  is  located  approximately  in 
the  middle  of  the  data  window.  Because  the  correlation  is 
modulo  P,  the  second  component  occurs  at  (P/2  +  0.75P)  mod  P  =' 
P/4.  The  relative  delay  between  components  appears  as  P/4 
instead  of  3P/4.  Thus,  components  which  were  delayed  by  3P/4 
could  appear  to  be  delayed  by  P/4  and  the  reverse  could  also 
occur.  The  safe  approach,  used  in  analysis  of  data  from  the 
experiment,  is  to  measure  the  spread  under  all  possible  rotations 
of  the  data  window  modulo  P  and  select  the  minimum  delay  under 

all  rotations.  This  approach  yields  slightly  conservative  values 

✓ 

for  the  average  delay  spread. 
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